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The fluctuation eigenmode problem of the nematic topological disclination line with strength 
±1/2 is solved for the complete nematic tensor order parameter. The line tension concept of a 
defect line is assessed, the line tension is properly defined. Exact relaxation rates and thermal 
amplitudes of the fluctuations are determined. It is shown that within the simple string model of 
the defect line the amplitude of its thermal fluctuations is significantly underestimated due to the 
neglect of higher radial modes. The extent of universality of the results concerning other systems 
possessing line defects is discussed. 

PACS number(s): 61.30.Dk, 61.30.Jf, 05.40.-a 

The fluctuation problem of the nematic disclination line is one of the central yet untackled problems in the field 
of liquid crystals, despite the fact that topological disclination lines are the most distinct and regularly observed 
characteristic of the nematic phase [1-4] and are reflected even in its name. Lately they have been playing a role as 
experimental objects in studies related to cosmology [5-7]. Controlled experiments with isolated lines or pairs of lines 
can be performed to study their motion [8,9]. Thermal fluctuations of disclination lines can be observed and analyzed 
simply by polarization microscopy [10,11]. 

In this paper, we present the solution of the full tensorial fluctuation problem of a straight wedge disclination 
line with strength ±1/2, in terms of finding the exact eigenmodes. We focus on the soft modes that correspond to 
string fluctuations (Fig. 1, inset) and challenge the validity of the string model of the disclination line [12]. The 
main results of our study are not limited to nematics, as our description does not include the physical system-specific 
elastic anisotropy or hydrodynamic flow. In particular, they do not depend on the winding number of the defect and 
can be directly applied to disclination lines in smectics-C and vortex lines in superfluids. The fluctuations of the 
strength 1 nematic disclination line were studied in the director description [13,14], precluding the string fluctuations. 
Moreover, the topological ±1/2 disclinations cannot be described by the Frank director elasticity, whereas the integer 
disclinations are known to be unstable. 

The director free energy per unit length — line tension — of a straight disclination line with strength s is 

T = 7rs 2 Kln- 1 (1) 
ro 

where K is the Frank elastic constant, R the system size and ro a microscopic cutoff. Hence, in a model the disclination 
line can be considered as a simple string under tension [12, p. 178]. The energy cost of its (overdamped) fluctuation 
modes (Fig. 1, inset), which can be attributed to the increase in length, Al — Jdz(du/dz) 2 /2, is thus 



AF = -k u To J dz cos kz, (2) 

where u is the amplitude and k the wave vector of the mode. The string model of a line defect is general [15-17] and 
not restricted to nematic disclination lines. 
■ Let us set the scene for the complete description of the disclination line in terms of the nematic tensor order 
parameter Q. Cylindrical coordinates (r, (f>, z) with corresponding r-dependent orthonormal base vectors {e r , e^, e z } 
will be used. The unperturbed disclination line (the ground state, as referred to from now on) coincides with the z 
axis. The ground state is z-independent, hence, the z-dependence of the eigenmodes is sinusoidal. In one (elastic) 
constant approximation, the free energy density f(Q, VQ) reads 

/ = iATrQ 2 + ±5TrQ 3 + ±C(TrQ 2 ) 2 + ±£Tr(VQ • VQ) (3) 
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and is invariant upon a homogeneous rotation of the Q-tensor. Hence, the Q-eigensystem rotates as ip — ipQ + (s — 1)0 
with respect to the above base vectors when we encircle a defect of strength s located at the origin; ip is the angle 
between e r and the tensor axis corresponding to the director in the far-field, while ipo is the free parameter of the 
defect configuration, corresponding to the angle between the director axis at <fi = and the x axis. There is no 
dependence on <j) other than this rotation, i.e., the scalar invariants of Q (the degree of order and biaxiality) are 
(^-independent — a "generalized" cylindrical symmetry. 

Defining another r-dependent orthonormal triad {ei,e2,e z }, where ei — e r cos?/> + e^smip and e 2 = — e r smi(j + 
e^cosip, in the unperturbed configuration the Q-tensor eigensystem coincides with this triad everywhere. Defining 
further the orthonormal symmetric traceless base tensors [18]: T = (3e z <g> e z — \)/V6, Ti = (ei ® §1 — e 2 <8> e 2 )/-\/2, 
T_i = (ei <g> e 2 + e 2 <g> ei)/\/2, T 2 = (e 2 <g> §1 + ei <g> e z )/V2, T_ 2 = (e 2 ® e 2 + e 2 <g> e z )/V2, the ground state 
Q-tensor components are (/(-independent, while the components of an arbitrary perturbation can be translated in <f> 
not changing the energy. Hence, the (^-dependence of the eigenmode components is sinusoidal. 

In the ocean of fluctuations one should first look for the Goldstone modes, as they yield families of slowly relaxing and 
spatially extensive "soft" fluctuation modes that can be easily observed. "Massive" fluctuations are less interesting, 
as they are short-living (wlOO ns) and localized («10 nm). Besides the homogeneous rotations of the Q-tensor, which 
in a deformed nematic are soft only in the one constant approximation, in a deformed system there exist additional 
nontrivial Goldstone modes, in our case the interest is in those corresponding to the displacement of the disclination 
line. Since translational degrees of freedom are absent in the analysis, the "displacement" is achieved by modifying 
(perturbing) the order parameter field. Modulating the displacement modes sinusoidally along z results in the string 
fluctuations of the disclination line (Fig. 1, inset). 
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FIG. 1. Radial eigenvalue A r (in units of 1/r) and the norm squared of the radial eigenfunctions (arbitrary units) for two 



radii of confinement. For R = 200, = 1.4 • 10~ 5 , = 4.5 ■ 10~ 4 , X^' = 1.4 • 10~ a ; note that X^'/X). 1 ' « 30! Inset: the 
string fluctuation of the disclination line. 

Expressing the Q-tensor as Q(r, t) = qi(r, t)Tj(r), the elastic (gradient) part of (3) is 

r = \ {(^) 2 + (^) 2 + 



(3) 



v (2) /\(1) 



(4) 



1 

77 



while the homogeneous (gradient-independent) part of (3) is a polynomial in qi. We introduce dimcnsionless quantities 
r <— r/£, t <— t/r, (A,B,C) <— (A, B,C)£, 2 /L, with the correlation length of the degree of order £, typically a few 
nm, and the characteristic time r = ^i^ 2 /i, typically tens of ns, where [ii is the bare rotational viscosity, i.e., 
71 = 95 2 /ii/2; 71 is the usual rotational viscosity and S the degree of order. Neglecting the hydrodynamic flow, the 
order parameter dynamics is governed by the time-dependent Ginzburg-Landau equation 



df 



— , 



(5) 
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where the time derivative vanishes in the ground state. 

Owing to the generalized cylindrical symmetry, the ground state involves only Qo(r) = ao(r)To + ai(r)Ti, as 
opposed to perturbations: AQ(r, t) — Xi{r,t)Ti. The fluctuation eigenmodes Xi, satisfying ±i = —Xxi, are sought by 
the ansatze 



%i = Ri, m (r) $>i(m(j)) sin(fcz) exp(— At), i — 0, ±1, 
Xi — Ri t n{r) <&i(n(j)) sm{kz) exp(— At), i = ±2, 



(6) 
(7) 



where $i = cos for i > and $i = sin for i < (the global angular phase and the z phase are arbitrary), m is an 
integer, whereas n = 1/2,3/2,5/2, ... is a half-integer due to the continuity and differentiability requirements (spinor 
symmetry of the base tensors)! In the one constant approximation, the sets of components (6) and (7) are not coupled, 
i.e., in-plane and out-of-planc fluctuations are independent. Furthermore, the z-dependence is fully decoupled, i.e., 
the eigenfunctions Ri (r) do not depend on k and the cross-section structure of the disclination line is not affected by 
the z modulation. 

We focus on the in-plane fluctuations. Putting the ansatz (6) into Eq. (5), the eigensystem for the radial functions 
Ri,m( r ) remains, with X = X r + k 2 , where A r is the radial eigenvalue: 



V Ro,m +[X r - go(r) - — ) i?0,m + 001 M Rl, m = 0, 



V 2 i? M „+ X r -gi{r) 



m 2 +4s 2 
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gai(r) i? ,m = 0, 



>-■' , * i \ m 2 +As 2 \ 4sm _ 

V R-l, m + [X r ~ .9-1 (r) - R-l,m o~~ R hm = 0; 



(8) 
(9) 
(10) 



the g's are quadratic polynomials of the ground state components. Note that the operator (8)-(10) is self-adjoint. 
Also note that defects with strengths s and — s are formally equivalent, i.e., changing the sign of the defect and 
redefining T/_ 1; _2} — > — T/_i _2} conserves the sign of s in the equations. 

For an infinite system, one can determine the homogeneous displacement mode directly by construction. Putting 
Qo(r — u) = Qo(r) + AQ(r), the perturbation AQ corresponding to the displacement u is AQ(r) = u • <9Q /<9r, 
which reads 



a;o(r)To + xi(r)Ti + z_i(r)T_i = -u • ( e r -^T + e r -^Ti + e 2syT_i 



(11) 



Thus, Xi in (11) are the eigenfunctions corresponding to A r = and m = 1. The lowest family of string modes is 
generated by adding the ^-dependence. Hence, A = fc 2 , or A = Lk 2 j \i\ = Kk 2 /ji in physical units (K = 9S 2 L/2), 
which is thus the exact result for the relaxation rate of the string mode, not obtainable in the director picture. 

The energy cost of an arbitrary fluctuation eigenmode is given by AF(x) = A JdVx 2 /2. For the homo geneous 
displacement mode (11), integrating over 0, this results in 



AF — -k 2 "- 2 



cos kz. 



(12) 



Comparing (12) (multiplied by Lt; to revert to physical units) with the form of the energy cost (2) of the simple string 
fluctuation, a line tension can be defined for the disclination line, as indicated by the underbrace. If the string model 
was accurate, the line tension defined in (12) would be the actual free energy (3) of the unperturbed disclination line 
per unit length. One verifies that it is exactly the elastic free energy (4) per unit length, whereas the homogeneous 
contributions are absent. Only the elastic terms contribute to the line tension! This is a general statement — any 
other homogeneous contributions, should they exist, e.g., electric, do not enter the line tension. It must be stressed 
that within the isotropic order parameter elasticity this finding is universal. Physically it means that the string 
fluctuation conserves the volume of the defect line, as opposed to the fluctuation of a real string. Thus, far from the 
defect core, where the homogeneous free energy contributions vanish, the string model is exact. Moreover, using the 
bulk values of a and a\, the line tension (1) is recovered exactly. Due to the director distortion, however, this regime 
is not approached exponentially but by a power-law. 

Let us now focus on the m = 1 fluctuations with A r > 0. Unlike R-i, the functions Rq and R\ are "localized" 
(Fig. 2), since go, gi — * const. > (of order 1), while g_i — > — 4s 2 /r 2 . Asimptotically, R-i, m (r) behaves as a 
combination of Bessel functions J m {V\-r) and Y m (y/X^r), whereas -R{o,i} nrs t decay as exp(— r)/y/r, followed by a 
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power-law asymptotics (Fig. 2, inset): R 



{o,i}.m ^ 4msfi_i >m /r 2 . In solving (8)-(10), we confine the system at a 

radius R with the restriction Ri >m ,n{R) = in order to get a discrete spectrum \i n \ n = 1, 2, 3, ... One might argue 
that this boundary condition is rather unphysical. In a real sample, however, one never deals with an ideally isolated 
disclination line — the observed line is surrounded by other defects and irregularities. The boundary condition should 
therefore be viewed as due to an effective confinement. In the limit R — > oo, physical observables should be only weakly 
dependent of R anyhow, e.g., like the free energy (1). Furthermore, if one uses the boundary condition R'^R) = 0, 
the lowest eigenmode is of the growing type, A r < 0, where A r oc 1/R 2 , reflecting the instability of the defect towards 
the escape from the system. This finite-size effect is often manifest in experimental situations. 
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FIG. 2. Three lowest radial eigenfunctions (not normalized) for m = 1. Functions Ro,i, n and Ri,i, n overlap on the scale of 
the big graph; the detail for n = 3 is shown in the inset — note the power-law tail and the zeros, which coincide with the zeros 
of R-i. 



The boundary conditions at r = are obtained by finding the analytic behavior of the ground state and the 
perturbations for small r: a' — > 0, a\ — > const, x r' 2s ' , Ra, rn — > r m , and i?±i, m — > ±a 1 r\ m ~ 2s \ +a 2 rl m+2s L The 
eigensystem (8)-(f0) is discrctized and efficiently solved by a multidimensional Newton relaxation method [19, p. 
588], Fig. 2. In the relevant limit R ^> £, the radial dependence of the lowest modes is identical for r <C R and 
given by (11), except for the normalization, i.e., each of these modes contributes to the displacement u of the central 

part! Still in the limit R » £, the difference \J Ar™ +1 ^ — \J A^ approaches ir/R when going to higher modes, i.e., 

Ar™^ — > (const. + mr/R) 2 , Fig. f. Similarly, for n ^> 1, the eigenvalue of the nth mode scales as A[™' oc f / R 2 with the 
system size R. 

We have shown that within the simple string model the degrees of freedom of the disclination are cut to the lowest 
m = 1 mode only. Thus, quite universally, the contributions of the higher modes represent a renormalization of the 
model. Thermal fluctuation amplitude at r — of the kth. Fourier component (z direction) is 

(AQ?(r = 0,m = 1, k)) = 2 £<c£,»>fl?i ,„(0), (13) 

n 

with the radial functions normalized, (c 2 n ) = IzbT '/ 'L£(k 2 + Ar™' 1 ) given by equipartion, and the factor 2 coming from 
the twofold degeneracy (cos</> and sin^ in the angular part). The displacement (u|) of the line is obtained from (13) 
by means of (11). If only the lowest radial mode is taken into account in (13), which is equivalent to using the non- 
renormalized string model, the thermal amplitude is obviously underestimated, Fig. 3, or conversely, the extracted 
value of K is too low. The error increases rather slowly (logarithmically) with R; for experimentally relevant scales 
it is of the order of 100%. The failure of the string model is best demonstrated in the limit R — > oo: for a fixed and 
nonzero k, (u|) — ► if only the lowest radial mode is taken into account, since the norm (squared) of the mode is 
logarithmically diverging when R — > oo. Alternatively, the same can be seen by noting the logarithmic divergence of 
the line tension (1). 
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FIG. 3. rms displacement of the disclination line (the length unit is £) vs. the number of radial modes summed in the regimes 
of a) weak and b) severe confinement; n = 1 corresponds to the non-renormalized string model. The maximum size of the 
computational domain is R = 580, the values beyond this range (R = 1000, R = 10000) are extrapolated. 



Thus, to simulate or analyze a real observation of the disclination line fluctuations, e.g., by polarizing microscopy, 
for a given k one has to sum up a proper number of radial modes possessing a considerable thermal amplitude; the 
maximum number is set by the resolution of the instrument — the summation should be stopped at the latest when 
the first zero of the radial functions becomes comparable to the resolution. As in a real situation 1/k is large compared 
to the resolution, the series is always truncated earlier by the diminishing thermal amplitude and hence there is no 
cutoff ambiguity. The two characteristic regimes are illustrated in Fig. 3: a), if ir/k <C R, many radial modes have to 
be summed up and the fluctuation amplitude is essentially independent of the system size R, whereas b), for ir/k > R 
(severe confinement) a few modes suffice and the fluctuation amplitude is suppressed by decreasing R. 

Of course, there exist other soft fluctuation modes, i.e., those with m ^ 1, which have not been focused on. 
Their radial eigenfunctions have zeros at r = 0, and hence, these modes are not important for the disclination line 
fluctuations, neither are they particularly characteristic for the defect structure. The m — soft mode, which involves 
only the component X-i (Eq. (10) is decoupled in this case), is a special case, corresponding to the Goldstone rotation 
of Q around the z axis. Naturally, it is present regardless of the configuration; the defect structure brings about 
merely a decay of the eigenfunction at the uniaxial core. 

Our results are valid in the approximation where the free energy is invariant to separate rotations of space and the 
order parameter (nematics and smectics-C in the one constant approximation, supcrfluids). More precisely, we require 
that the order parameter gradients in the z direction enter the free energy in the form of square terms only, Eq. (4) . If 
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this symmetry is broken, the gradients along z arc in general coupled to those in the xy plane. In a related manner, 
superconductors and their analogues, smectics-A (SmA) [20], besides the AF-model (nematics, superfluids) degrees 
of freedom '5 = |\I/|expi(/> (wave function of the superconducting state or the amplitude and phase of the smectic 
modulation) possess an additional vector order parameter (vector potential A or the c-director, respectively). Due to 
the coupling between V^P and A described by the free energy density term | (-iW — eA) ^| 2 /2m, gradients in the z 
direction induce a change in A (analogously for SmA). As a consequence, in the string fluctuations the cross-sectional 
structure of the line defect (vortex or dislocation) is not merely displaced but also deformed, and, hence, the actual 
line tension is not simply the free energy per unit length (not even asymptotically). In this case one can define an 
effective line tension and must still renormalize the string model with the contribution of higher radial modes. 

This work was supported by the Slovenian Office of Science (Program Pl-0099) and US-Slovene NSF Joint Found 
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